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ABSTRACT 

There is no universally acknowledged criterion to distinguish brown dwarfs from planets. 
Numerous studies have used or suggested a definition based on an object's mass, taking the 
~13- Jupiter mass (Mj) limit for the ignition of deuterium. Here, we investigate various 
deuterium-burning masses for a range of models. We find that, while ISMj is generally a 
reasonable rule of thumb, the deuterium fusion mass depends on the helium abundance, the 
initial deuterium abundance, the metallicity of the model, and on what fraction of an object's 
initial deuterium abundance must combust in order for the object to qualify as having burned 
deuterium. Even though, for most proto-brown dwarf conditions, 50% of the initial deuterium 
will burn if the object's mass is ~(13.0 ± 0.8)Mj, the full range of possibilities is significantly 
broader. For models ranging from zero-metallicity to more than three times solar metallicity, 
the deuterium burning mass ranges from ^11.0 Mj (for 3-times solar metallicity, 10% of initial 
deuterium burned) to ^16.3 Mj (for zero metallicity, 90% of initial deuterium burned). 

Subject headings: radiative transfer - stars: low-mass, brown dwarfs - stars: evolution 



1. INTRODUCTION 

The year 1995 heralded both t he first u n ambigu- 
ous d etection of a brown dwarf (Oppcnheimcr et al.l 
119951 ) and the first unambig uous detections of plan- 
ets beyond our solar s ystem (jMavor fc Quelodll995l : 
iMarcv &: ButleHll996D . Many of the first substellar 
objects detected were either clearly brown dwarfs 
(very massive, not in a close orbit around a main- 
sequence star) or clearly planets (lower mass, orbit- 
ing close to their stars). However, as the number of 
discoveries of substeller objects grew to the dozens 
and then hundreds, there increasingly appeared to be 
overlap in the apparent mass distributions of brown 
dwarfs and planets. This highlighted the need to 
clarify the taxonomy. Moreover, the various defini- 
tions have been strained by the recent discoveries of 
objects such as CoRoT-3b, a ~22- Jupiter-mass (Mj) 
object in a close ( 0.057-AU) orbit around its star 
(jDeleuil et al.ll2008 ) and the directly imaged objects 
of masses ~ 5-15 Mj at t ens of AU from HR 8799 and 
Fomalhaut (jMarois et al . 2008; Kalas_et_al, 2008). 

A commonly used way to classify objects that are 
~10-15 times the mass of Jupiter is by whether they 
fuse deuterium (D) in their deep interiors. This cri- 
terion was adopted in 2002 by the Working Group on 
Extras olar Planets of the International Astronomical 
Union (|Boss et al.ll2007| ): 

1 ) Objects with true masses below the 

limiting mass for thermonuclear fusion of 

deuterium (currently calculated to be 13 

Jupiter masses for objects of solar metal- i indeed, note that the lAU did not adopt this criterion. 

See, e.g .. _ 

I dsp@astro.princeton.edu, burrows@astro.princeton.edu, milsomP]l)tp:^pa gbrialirtatlea dty|edu/'^basri/defineplanet/Mercury.htm| 



licity) that orbit stars or stellar rem- 
nants are "planets" (no matter how they 
formed). The minimum mass/size re- 
quired for an extrasolar object to be con- 
sidered a planet should be the same as 
that used in our Solar System. 

2) Substellar objects with true masses 
above the limiting mass for thermonu- 
clear fusion of deuterium are "brown 
dwarfs, " no matter how they formed nor 
where they are located. 

3) Free-floating objects in young star 
clusters with masses below the limiting 
mass for thermonuclear fusion of deu- 
terium are not "planets, " but are "sub- 
brown dwarfs" (or whatever name is 
most appropriate). 

Although we (and others) do not necessarily en- 
dorse the "deuterium-burning edge" as the most use- 
ful delineation bet ween p lanets and brown dwarf^H 
fsee iBurrows et all l200lt iChabrier et"all 120051 and 
iBakos et al.ll2010f r it is a commonly used criterion 
and warrants further exploration. 

Models of brown dwarfs and giant planets have 
been calculated since before the first such objects 
were discovered ("Kumai^'igGSl: iZapolskv fc Salpetei 
1969 ; IGrossman & Grabos kT 119731: IBurrows et al ' 
1993t ISaumon et all il994t iBurrows et al.l I1995D 



Such computations depended heavily on the equa- 
tion of state (EOS), and were therefore aided by 



2 



the p ubUcation of a new EOS for H2-He mix- 
tures (jSaumon et all I1995D (the so-called "SCvH" 
EOS). In addition, models of the emergent ra- 
diation from, and of the temporal evolution of, 
brown dwarfs were improved by t he application 
of nongray rad iative transfer t heory (Burrows et alJ 
1997. 1998; B urrows fc Shard ir999: Burrows 1999; 

20011 I2 003r Tiaraffe et al. 2002, 
2003; Chabrier et al. 2000b, a; 



Burrows et al 
2003; All ard et al.l 
Chabrier fc Baraffd 



200d ISharp fc Burrowsl 12007^ 



Barmani" 2008f) and sop histic ated chemical rnod 



2002 : 



2007 



els (iFeglev &: Lodder^[T9 94: Lo dders &: Feelevl 
Burrows & Sharp" '199f* iSharp fc Burrowsl 
Allard & Hauschildt 1991). 

The theoretical study of brown dwarfs and mas- 
sive planets is now a maturing field that is of particu- 
lar current interest, given the increasing pace of dis- 
covery of objects in this mass range. Previous model 
calculations have suggested that deuterium burning 
turns on near a mass of 13 Mj ( Burrows et al.lll997l 
20011) or at 0.012 Mq (-12.5 Mj. iChabrier et al l 



20051) R An early such calcula tion was performed by 



Grossman fc Graboskd ()1973[ ). who claimed the ex- 



istence of a deuterium main sequ e nce ne ar O.O12M0, 
although ' Dantona fc Mazzitellil ()1985[ ) found, us- 
ing a lower deuterium abundance, that there is 
no deuterium main sequen ce. The models of 
iGrossman fc Graboskd (|1973l l used a deuterium frac- 
tion that was about 10 times the Galactic value, 
and the atmosphere treatment was less sophisticated 
than the current state-of-the-art. Although previ- 
ous calculations, and community prejudices, have 
converged on a mass limit of ~13 Mj, and this 
mass range does provide a reasonable estimate for 
the mass at which significant deuterium burning oc- 
curs, it is worthwhile to clarify the precise range 
of deuterium-burning-mas s predictions under a vari - 
ety of model assumptions. iSaumon fc MarlevI (|2008f ) 
have contributed to this endeavor. Specifically, defin- 
ing the minimum deuterium-burning mass as objects 
that burn 90% of their initial deuterium in 10 Gyr, 
they find 13.1 and 12.4 AIj for the cloudless and 
cloudy cases, respectively. This shows that how one 
computes the radiative boundary condition (with or 
without clouds, and what type of clouds), matters 
significantly. 

In this paper, we systematically explore how deu- 
terium burning depends on a several physical prop- 
erties of objects. Our goal in this paper is to point 
out that the division between objects that burn deu- 
terium and those that do not is not strictly 13 Mj. 
One important aspect of this is that deuterium burn- 
ing does not "turn on" suddenly at a particular mass. 

Calculations such as this implicitly assume that objects 
start with large initial entropies (which leads, eventually, to 
high central temperatures) , as is expected for bodies t hat form 
from the collapse of a cloud of molecular gas ( Marl ev et al.l 
I2007i ) . Note that if the initial entropy is low, objects that are 
significantly more massive could, in principle, avoid burning 
any significant amount of deuterium, though formation sce- 
narios such as this are probably not possible. 



Even with just a single set of model parameters, 
the mass at which 10% of deuterium burns is dif- 
ferent from the mass at which 50% or 90% burns. 
This ambiguity is only enhanced when one exam- 
ines the range of possible or likely values of the 
variables under consideration here (such as helium 
abundance, initial deuterium abundance, and metal- 
licity). We consider various (cloud- free) models and 
nuclear burning criteria; in ^ we describe the mod- 
els under consideration. We vary the helium fraction 
as a proxy for the effect of varied concentrations of 
elements heavier than hydrogen ( H2.2\\ : we vary the 
initial fraction of deuterium ( ij2.3p : and we vary the 
metallicity in the context of its influence on atmo- 
spheric opacity and cooling f iJ2.4l) . In we de- 
scribe the influence that each aforementioned knob 
has on deuterium burning, and in 52J we calculate 
derivatives of the D-burning mass edge with respect 
to each quantity (Y, j/d.i, Z). In ij5l we conclude 
by discussing the implications of our results for the 
distinction between brown dwarfs and giant planets. 

2. MODELS 

We consider a range of brown dwarf/planet mod- 
els. If the objects are massive enough, their deep 
interiors have conditions under which non-negligible 
burning of deuterium occurs through the reaction 
p + d — >■ 7 -|- '^He (by far the most important re- 
action involving deuterium at the temperatures and 
densities in the cores of these objects). The rate of 
deuterium burning depends on the concentration of 
deuterium, th e mass density, and the temperature. 
iStahleil ()1988f ) presents a power-law approximation 
to this rate in the range T ~ 10*^ K: oc pyj^T^^-^, 
where ed is the specific deuterium-burning rate, p 
and T are central values of density and tempera- 
ture, and ud is the deuterium mixing ratio. Al- 
though this expression does not, of course, describe 
the actual deuterium-burning rate used in our code, 
it may serve as a rough guide. We vary the helium 
mass fraction {Y, which influences the central den- 
sity), we vary the initial deuterium mixing fraction 
(2/D,i), and we vary the metallicity of the atmosphere 
{Z, which strongly affects the opacity). Our models 
are described in ij2. 1112.41 below and summarized in 
Table [T] Note that each row of this table actually 
corresponds to a family of models with a dense spac- 
ing of masses, ranging from 10 to 20 times Jupiter's 
mass. 



2.1. Modeling Technique 

The basic modeling strat egy is the one de- 
scribed in lBurrows et al.l H1997D . Specifically, we self- 
consistently calculate the evolution of these objects 
by matching radiative-convective atmosphere mod- 
els onto fully convective interiors. Objects are as- 
sumed to be born high on the Hayashi track, with a 
given high interior entropy and a flat entropy profile. 
This is a reasonable assumption because in Nature 
these objects become fully convective and lose their 
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initial conditions (i.e., the entropy decreases rapidly 
from its initial high value) within 10^~^ years af- 
ter formation. This is far shorter than the duration 
or age of deuterium burning, which can last lO''"^ 
years. Interior structure is calculated with the SCvH 
EOS (jSaumon et al.|[T995l) . In the atmosphere, con- 
vecti ve heat flux is treated with mixing length the- 
ory ()Mihala£lll978D . takin g the mixing le ngth to be 
the pressure scale height. iBurrows et al. l (1989) in- 
vestigated how other choices of mixing length affect 
the thermal evolution of brown dwarfs and low-mass 
stars and found that, for brown dwarfs, unlike for 
stars on the main sequence, this dependence is very 
weak. 

At temperatures found in the interiors of these 
objects, the rate of proton-deuteron fusion is re- 
ported in the literature as the standard extrapolation 
from the rates at higher tem peratures . We use the 
rate reported by iCaughlaiT fc Fowlei] (|1988l). w hich 
is the same as that found in lHarris et al.l (|l983l ). In 
1999, the Nuclear Astrophysics Compilation of RE- 
action rates (NACRE) released a new com pilation 
of thermonuclear rates (lAngulo et al.lll999f). which 
was u pdated several vears later (|Descouvemont et al.l 
l2004f ). S till, the recent rele ase of the Reaclib 
Database (jCvburt et al.l 120101 ) by the Joint In- 
stitute for Nuclear Astrop hysics (JINA) uses the 
iCaughlan fc Fowled (|1988l ) rate. A good dis- 
cussion of nucle a r fusi on rate s can be found in 
lAdelberger et al.l (l2010l). The IC aughlan fc Fowleii 
rate and that of iDescouvemont e t al. (2004) differ 
by '--^20% at T ^ 10^ K, and start to converge at 
lower temperatures. The interiors of these substellar- 
mass objects reach temperatures of ^^0.5-1x10^ K, 
where the Gamow peak energies are ~0.7-l keV. 
At these temperatures and energies, the discrep- 
ancy between the older and the newer estimates 
of the rate is less than 10%, although extrapola- 
tions in this regime are generically problematic. To 
our knowledge, there is no newer information on 
how properly to make this extrapolation to temper- 
atures and Gamow peak energies of relevance to the 
question at hand0 Screening corrections for both 
ions and electrons are emploved (ISaumon et "arill996t 
IDewitt et al1ll996l:ISahriing fc Chabriedll998D . 

The thermal evolution depends on the cooling 
rate as a function of o bject size and entropy. In order 
to calculate this, as in IBurrows et al.l (jl997D . we pre- 
calculate a large grid of one dimensional, non-gray, 
radiative-convective atmosphere models with differ- 
ent effective temperatures (Tcs) and surface grav- 
ities (g). The effective temperature specifies the 
net radiative-convective flux at every vertical level 
{aT^ff, where a is the Stefan-Boltzmann constant). 
Each model in the grid is fully convective at its base. 

^ The ambiguity in the proper prescription for deuterium 
fusion, of course, translates into some ambiguity in the 
deuterium-burning mass beyond what is suggested by the pa- 
rameter study described in this paper; this effect, however, is 
small. 



The grid of atmosphere calculations, therefore, can 
be thought of as defining a function S'[rcff, g], where 
S is the entropy. What is important for the thermal 
evolution is the net cooling rate (essentially Tcff). 
Therefore, we invert this function to obtain Tcs [S, g] . 
This function depends on the atmospheric opacities 
and is what is meant by the "surface boundary con- 
dition." At every step of the thermal evolution, the 
cooling rate is calculated by interpolating in the grid. 
In this way, the fully convective interior matches 
smoothly onto the appropriate atmospheric bound- 
ary condition. 

2.2. Helium Fraction 

For an ideal gas, density depends on molecular 
weight, while for degenerate matter, density is in- 
versely related to the ratio of electrons to baryons. 
The abundance of heavier-than-hydrogen elements 
affects the density for both reasons. We vary the 
helium content of our models from Y — 0.22 to 
Y — 0.32, with more helium-rich models represent- 
ing objects that are richer in both helium and metals. 
The corresponding model names (in Table [T] and in 
the figures) are He22 - He32. Model He25 has he- 
lium mass-fraction 0.25, and is also referred to as 
model D2 (described in H2.'3\\ . These models have a 
radiative boundary condition (i.e. at mospheric opac- 
ities) calculated using the method of IBurrows et al.l 
([1997 ') ■ with solar- metallicity opacities and initial 
deuterium number fraction of yu,i ~ 2 x 10^^. 

We note that not all helium abundances are 
equally likely. Historically, quite low values of Y 
have been inferred for Jupiter and S aturn. Using 
two d ifferent instruments on Voyager, IGautier et al.l 
(|1981l ) estimated that the helium mass-fraction of 
Jupiter was 0.1 9 ± 0.05 or 0.21 ± 0.06 , respectively. 
More recently, Ivon Zahn et al.l (|1998D used Galileo 
data to find the helium mass-fraction of Jupiter to 
be 0.234 ± 0.005. In a re analysis of Voyager data, 
IConrath fc Gautieil (j2000D found the abundance in 
Saturn to be in the range 0.18-0.25. One must 
keep in mind that, because of differentiation, atmo- 
spheric values for these objects might be lower than 
bulk abundances. Furthermore, the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) indicates that 
a fairly stringent lower bound on th e cosmic pri- 
mord ial helium mass fraction is 0.248 (jSpergel et al.l 
120071 ). Higher values than this can easily be ex- 
plained in some objects from Galactic enrichment, 
but lower values might be difficult to explain. Never- 
theless, we examine the effect of very low values oiY, 
down to 0.22, so as to explore the dependence of the 
deuterium-burning edge on Y values that have been 
invoked in the last several decades, but the reader 
should note that such an extreme value might not 
obtain in actual substellar objects. 

2.3. Initial Deuterium Fraction 

Big-bang nucleosynthesis calculations, in con- 
junction with other observations (CMB, high-z 
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quasars, etc.) suggest that the primordial D:H ratio 
is (2.8±0.2) X 10-5 (IPettini et aL|[2008l ). Our galaxy, 
however, is somewhat depleted in deuterium, relative 
to this value. Recent observations indicate that the 
mean mixing ratio of deuterium in the interstellar 
medium is roughly (2.0 ± 0.1) x 10"^, though the 
scatter (among different lines o f sight) is a factor of 
several ([Prodanovic et aLll2010t) . In our models, we 
take initial deuterium fractions of 10"^, 2 x lO'^, and 
4 X 10-5 (mixing ratios) The corresponding model 
names are Dl, D2, and D4. These models also have a 
radiative boundary condition (i.e. at mospheric opac- 
ities) calculated using the method of iBurrows et al.l 
([1997) , with solar- metallicity opacities and a helium 
mass-fraction of 0.25. 

2.4. Metallicity 

For a non-irradiated object, the metallicity, i.e. 
the abundance of elements heavier than helium, 
is the main determinant of the atmospheric chem- 
istry and, hence, the opacityQ Among objects 
with thick H2-He-dominated atmospheres, this abun- 
dance can range from very low, to more than 30 
times the solar abundance, as is expected for ob- 
jects such as Uranus, Neptune, and similar-mass 
exop lanets (Guillot & Gautier '2007: Spiegel et alj 
201131: Lewis et al. 2010; Ncttelmann et al. 2013 
Madhusudhan fc SeageHiMlOf ). StiU, objects as mas- 
sive 10 Mj might generally be expected to have 
metallicities not much more than 0.5-1 dex more 
than solar. 

In our models, we take metallicity values rang- 
ing from zero up to a half dex more than so- 
lar, and we calculate th e opacity corresponding 
to equilibrium chemist ry (jBurrows fc Shard IT999[ 
ISharp fc Burrow^ 120071 ). Here, we take solar metal- 
licity to be defined by lAnders & : Grevesse (1989). 
The typical derived atmospheric molecular opacities 
and, hence, the net cooling rate resulting from these 
older abundances generally do not di ffer much from 
those calculated using the more recent lAsplund et al.l 
(|200i^ values, and the latter remain controversial. 
Furthermore, we point out that using the older ele- 
mental abundances is no obstacle to our basic goal, 
which is to explore the broad consequences of differ- 
ent metallicities on thermal evolution. Since other 
stars have different elemental ratios than the Sun, 
the degree of uncertainty in how our mix of el- 
ements matches that of typical stars depends at 
least as much on variability in typical stellar val- 
ues as on the estimate of solar abundances. In 
the zero-metallicity models, atmospheric opacity re- 
sults entirely from coUisionally-induced absorption 
due to the collisional ly-induced dipoles of II2 and He 
(jBorvsow et al.l [1989 1. Though the bottom end of 
this range (0 and 0.01 times solar) includes lower 

V alues below this range are not rare l lProdanovic et all 
[2OT0l 'l. 

^ Metals affect interior opacity too, but at any metallicity 
brown dwarf interiors arc fully convective. 



metallicities than might not be expected in most 
brown dwarfs, it is interesting to consider the full 
range of possibilities. 

The models in which we vary metallicity fall into 
two sequences, ZO.Ol, ZO.l, Z0.3, Zl.O, and T0.3, 
Tl, T3. Models in the form er sequence are cal- 
culate d using the opacities of lAllard fc Hauschildtl 
( 1 99-"h to derive the atmospheric boundary con- 
ditions (D. Saumon, private communication) that 
drive the evolution. These models have a helium 
mass-fraction of 0.25 and initial deuterium fraction 
of 2 X 10^^. The metallicity of these models (as a 
fraction of solar) is the number following the Z in 
the model name. Those in the latter sequence have 
a boundary co r idition calculated using COOLTLUSTLY 
(IHubenvl 119881: IHubenv fc Lanzlfl99 a IHubenv et all 
120031: IBurrows et al.l 1200 6: Sharp & Burrowsl 12001 7 
with the same helium and initial deuterium fractions 
(|Heng fc Burrowsll20T0 in prep. I ). Among these, the 
metallicities of TO. 3 and T3 are not exactly 0.3 and 
3, but instead 10"^/^ and 10^/^ times solar metal- 
licity. Finally, model ZO has zero metallicity and is 
grouped by label with the former sequence, but actu- 
ally has a rad i ative boundary condition taken from 
ISaumon et all ()1994D . 

3. RESULTS 

By considering a range of model parameters {Y, 
yu.'n Z) and atmospheric radiative boundary con- 
ditions, and by examining evolutionary trajectories 
for a range of masses, at ages from 1 Myr to 10 Gyr, 
we have produced a broad range of model calcula- 
tions. Here, we show how properties such as ra- 
dius, effective temperature, and the deuterium burn- 
ing rate and power vary with time as a function of 
mass (^ ii3.ip . Furthermore, we examine how these 
evolutionary tracks vary with the other parameters 
at fixed mass f t|3.2p . and we calculate the variation 
of the deuterium-burning mass limit with the model 
parameters listed above fi i|3.3|) . 

3.1. Evolution with Mass 

The basic character of brown dwarf cooling tra- 
jectories is well known from a number of calculations 
in the last several decades, but it is worthwhile to ex- 
amine the specific behavior of our models. We begin 
with an illustrative set of evolutionary trajectories 
for a representative family of models (He25, or D2). 
Figure [T] shows the evolution of four quantities that 
our models track. 

The top panels show quantities related to model 
cooling. The top left panel shows the evolution of 
radius with time, and the top right panel shows the 
evolution of effective temperature with time. It is 
instructive to examine these two panels together. 
Higher mass objects begin larger but have higher ef- 
fective temperatures and therefore cool more rapidly. 
As a result of their faster cooling, the radii of more 
massive objects shrink more rapidly, eventually (af- 
ter roughly a gigayear) overtaking lower mass ob- 
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jects in radius. By the end of the calculation (after 
10 Gyr), all objects are very nearly Jupiter's radius, 
but more massive objects are slightly smaller. 

The lower panels of Fig.[l]are related to the fusion 
of deuterium. During the initial Hayashi-track stages 
of evolution, the luminosity is dominated by Kelvin- 
Helmholtz contraction. After a few million years, 
the central part of the object reaches temperatures 
and densities at which the deuterium fusion rate is 
non- negligible. The bottom left panel shows the ra- 
tio of the power from nuclear (deuterium) burning 
to the total luminosity. For the most massive ob- 
jects, this ratio reaches fairly large values (well over 
80%). However, the main sequence is where this ra- 
tio is essentially 100% (i.e., where Kelvin- Helmholtz 
contraction has ceased), and so it is clear that none 
of these models goes through a deuterium main se- 
quence stage. The bottom right panel shows the evo- 
lution of the fraction of the initial deuterium (in this 
case, j/D,i = 2 X 10~^) that has burned. Higher mass 
objects burn a larger fraction of their initial deu- 
terium, and do so faster than lower mass objects. 
Regardless of mass, very little deuterium burns af- 
ter a few hundred million years. Note, in compar- 
ing the top row to the bottom row, that the 15-Mj 
and 20-Mj models, which fuse the greatest amount 
of deuterium (among the models displayed), cool 
and shrink very quickly after the deuterium-burning 
phase. 

3.2. Model-Dependence of D-Burning Evolution 

In Fig. [21 we examine the deuterium-burning evo- 
lutionary profiles of 13- Mj objects from a variety of 
models. Each row shows analogous plots to the bot- 
tom row of Fig. [l] only instead of different curves 
representing different mass objects, they represent 
different model properties. The left column shows 
the evolution of the ratio of nuclear power to total 
luminosity, and the right column shows the fraction 
of the initial deuterium that burns. The three rows 
show, from top to bottom, models He22-He32 (vary- 
ing helium fraction), models D1-D4 (varying ini- 
tial deuterium fraction), and models ZO-Zl.O (vary- 
ing metallicity through its influence on atmospheric 
opacity). Each row also shows models He25 (which 
is also D2) and Zl.O, for comparison to "fiducial" 
models. 

First, it is important to note the difference at fixed 
model parameters of the radiative boundary condi- 
tion. The Z-sequence models have higher opacity, 
which leads to slowe r cooling, than the other mod- 
els (which have the iBurrows et al.l [1997 boundary 
condition). As a result, models Zl.O and He25/D2, 
which are otherwise identical, have markedly differ- 
ent evolutionary trajectories. The former reaches a 
peak nuclear-to-total power ratio of «80%, while the 
latter reaches only Ri50%. Furthermore, the former 
eventually burns 66% of its initial deuterium, while 
the latter burns only 17% of its. In fact, of the Z- 
sequence models, model Z0.3 (with 0.316 x solar el- 



emental abundance) more nearly approximates the 
He25/D2 evolution than does Zl.O. With its slower 
cooling, model Zl.O has a moderate nuclear-burning 
to total luminosity ratio until later times; the ra- 
tio shown in the left column stays above 5% until 
-5 Gyr. 

More generally, within each row of Fig. [21 the 
expected trends hold. The top row shows that 
more helium-rich models burn a greater fraction of 
their initial deuterium, which is unsurprising because 
these are the models with fewer electrons to support 
the mass via degeneracy pressure; thus, these mod- 
els have denser, hotter cores. The middle row shows 
that models with higher initial deuterium content 
burn a greater fraction of their initial deuterium. 
This is a somewhat subtle point. The deuterium- 
burning rate is roughly linear in suggesting that, 
other things being equal, the fractional depletion of 
deuterium with time should be constant irrespec- 
tive of initial deuterium abundance. However, other 
things are not equal. In particular, the central tem- 
perature of more deuterium-rich models is main- 
tained at a higher value for a longer period of time. 
As a result, the integrated fractional amount of deu- 
terium burned is greater for these models. The bot- 
tom row shows that more metal-rich models burn 
a greater fraction of their initial deuterium, which 
stands to reason because their higher opacities pro- 
duce a blanket effect, allowing them to maintain 
higher core temperatures for a longer time. 

3.3. Model- Dependence of D-Burning Limit 

We now quantify the deuterium-burning mass 
limit, and its dependence on the various model pa- 
rameters we have varied. Figure [3] displays the 
same basic information in two different ways. In 
the left column, the fraction of the initial deuterium 
that combusts within 10 Gyr is plotted versus ob- 
ject mass. Horizontal dashed lines are plotted at 
10%, 50%, and 90%, illustrating various (arbitrary) 
deuterium-burning cut-offs. The right column plots 
the deuterium-burning mass edge for each of these 
three criteria (i.e., the mass at the intersection of 
each curve in the left column with the corresponding 
dashed line) as a function of the three tunable model 
parameters: Y, j/D,i, and Z. From top to bottom, the 
rows show the same respective model sets as in Fig. [2] 
As in Fig. [21 models He25/D2 and Zl.O are shown 
in all three panels of the left column. The last three 
columns of Table [T] contain the same information as 
the right column of Fig. [3l and also include data 
for the COOLTLUSTY sequence of models (T0.3-T3). 
Note that the solar- metallicity COOL T LUSTY model 
(Tl) and the fiducial IBurrows etHI (|1997[ ) model 
(He25/D2) are quite similar to one another. 

There are four main lessons from the data in 
Fig. [3] and Table [B 

1. Within each sequence of models, the ex- 
pected trends hold. That is, greater helium 
abundance, greater deuterium abundance, and 
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higher metalhcity aUow a given amount of deu- 
terium to fuse at a lower object mass. 

2. Depending on what is meant by "deuterium 
burning" (i.e., how much of the initial deu- 
terium must burn to qualify), the deuterium- 
burning mass limit can vary quite a bit. As 
the criterion goes from 10% to 50% and from 
50% to 90%, the required mass increases by 
~0.8-0.9 Mj. 

3. 0.012 Mq and 13 Mj is not far from the mass 
limit for most models. 

4. However, the full range of masses displayed in 
this table goes from 11.3 Mj (He32, 10%) to 
16.3 Mj (ZO, 90%). 

With the final point in mind, it is clear that claiming 
that an object does or does not "burn deuterium" im- 
plies a set of physical assumptions about the model 
and a criterion for deuterium burning. Any criterion 
is fine, so long as people know what it is; but implicit 
assumptions should be made explicit, and might not 
actually apply to some astrophysical objects. 

4. DISCUSSION 

A useful way to quantify the influence of varying 
the model parameters on the deuterium mass limit is 
with the derivative of the mass-limit with respect to 
each model parameter. Although the curves in the 
right column of Fig. [3] do not have constant slope, 
they are not very far from linear over the ranges 
shown. In Table [2l we present the average derivative 
with respect to each parameter. For Y and utia, 
we take linear derivatives; for j/D,i and Z ^ we take 
logarithmic derivatives. These derivatives are not 
terribly precise measures of the variation, but allow 
for a quick, crude estimate of the deuterium edge- 
mass at different parameter values. 

Note that, in the calculations discussed so far, 
the effect of higher metallicity has been reflected 
only in the atmospheric opacity. In a real object, 
however, higher metallicity throughout would also 
lead to an EOS change for the bulk, because high- 
Z elements have lower electron-to-baryon ratios. As 
a result, a higher metallicity object would have a 
higher central density than occurs in our models, 
and would therefore have a lower deuterium mass 
limit than we have shown, particularly for model 
T3. Though there is no published robust equation of 
state that properly includes heavy elements beyond 
helium, we can approximately correct for this by us- 
ing a h i gher helium fraction (this was also done by 
iGuillotI 120081 ). We, therefore, have post-processed 
our models in an attem pt to incorpor a te thi s ef- 
fect. Using values from lAsplund et all (|2009l ). we 
take Y = 0.2703 and Z = 0.0142 to be the he- 
lium and heavy-element mass-fractions correspond- 
ing to protosolar abundance. Motivated by this, we 
add 0.0142 in Y for every unit in metallicity (i.e., 
Y = 0.2703+0.0142 for solar, Y = 0.2703+2x0.0142 



for 2x solar, etc.)I3 We perform a cubic interpo- 
lation in Y and use the ^'-derivative for model se- 
quence T0.3-T3 to calculate the values in Table |3] 
corresponding to 0.5, 1, 2, and 3 times solar metal- 
licity. Values of the edge mass in this table then 
range from ~11.0 Mj (3x solar, 10% of initial deu- 
terium burned) to -13.9 Mj (0.5x solar, 90% of 
initial deuterium burned). 

Although the deuterium-burning mass limits in 
Tables [1] and [3] vary from 11.0 to 16.3 times Jupiter's 
mass, a narrower range of values is found if we re- 
strict our attention to a particular burning fraction 
criterion (50%, say) with realistic Y, j/d,!, and Z 
values. For helium mass-fractions between 0.25 and 
0.30, initial deuterium fractions between 10~^ and 
2 X 10~^, and metallicities between 0.5 and 2 times 
solar, the mass required to burn 50% of the initial 
deuterium is between ~12.2 Mj and ~13.7 Mj. Val- 
ues somewhat outside this range are not impossible, 
but values very far from this range are probably rare. 

5. CONCLUSIONS 

We have calculated a suite of models of substel- 
lar mass objects, encompassing a range of values 
of helium fraction, initial deuterium fraction, and 
metallicity. We have also used several calculations 
of the atmospheric boundary conditions. Although 
the rate of deuterium-burning is extremely sensitive 
to temperature, it is worthwhile bearing in mind that 
deuterium will burn at some rate at any nonzero 
temperature, so one must specify what is meant by 
"deuterium burning" (i.e., how much of the initial 
deuterium must burn to qualify). We have found 
that, depending on what criterion is used, reason- 
able values for the deuterium mass limit range from 
~11.4 to —14. 4 times Jupiter's mass; the former limit 
corresponding to 2x solar metallicity, 10% of ini- 
tial deuterium burned, and the latter limit to so- 
lar metallicity, Y — 0.25, 90% of initial deuterium 
burned. Extreme models (very low helium frac- 
tion, very low metallicity) could require even higher 
mass in order to burn a specifled percentage of their 
i nitial deuterium. T he canonical value of 13 Mj 
(jBurrows et al.|[l997l ) is a reasonable estimate of the 
mass at which signifcant deuterium burning begins, 
but this value is model-dependent. 

Still, one should keep in mind that the deu- 
terium cut is probably less relevant to an object's 
true taxonomic status than is its formation history. 
The downside of a "formation scenario" deflnition is 
that formation history is not easily observable. On 
the other hand, though mass might be observable, 
the parameters Y ^ yu.'n and Z might not be. Fur- 
thermore, some objects will be found significantly 
above 13 Mj that will be of clearly "planetary" ori- 
gin (e.g., CoRoT-3b). Whether to call these objects 
"brown dwarfs" or "deuterium-burning planets" re- 

^ It is not clear that the equivalent Ay for solar metal- 
licity should be exactly Z, but this seems a reasonable first 
approximation. 
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mains open for debate. Similarly, objects might be 
found that are significantly less massive than 13 Mj 
that appear not to have formed in a protoplanetary 
disk, but rather to be the low-mass end of the brown 
dwarf formation process. Though classified as "plan- 
ets" by standard terminology, these objects might be 
more taxonomically related to brown dwarfs than to 
planets. 

We emphasize that there is really no need at this 
time to rigidly distinguish between a giant planet 
and a brown dwarf on the basis of a single criterion. 
There is ambiguity in the provenance of these ob- 
jects, and this ambiguity might persist for a while. 
The use of a particular term ("planet" or "brown 
dwarf"), should be accompanied by the definition 
that is employed. Given that planets are thought 
to be objects in orbit around a star (or aroimd a 
brown dwarf), while brown dwarfs are thought to 



be the low-mass end of the star formation process, 
there is likely to be overlap in the mass range of 
these objects unless one adopts a rigid mass cut to 
distinguish them. Doing so, however, will surely lead 
to an overlap in formation scenarios (as in the case 
of CoRoT-3b). When classifying a newly-discovered 
substellar-mass object, one can use a variety of "rea- 
sonable" criteria, but the classifications should re- 
main tentative until a more thorough observational 
and theoretical understanding of substellar objects 
is achieved. 
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TABLE 1 
Models 



Model 


Z 


yo.i 




Y 


M (10%) 


M (50%) 


M (90%) 












(A/./) 


{Mj) 


{Mj) 


He22 


1 


2 X 10~ 


5 


0.22 


13.20 


14.08 


14.92 


He25 (D2)'" 


1 


2 X 10" 


5 


0.25 


12.70 


13.55 


14.35 


He28 


1 


2 X 10" 


5 


0.28 


12.06 


12.83 


13.62 


He30 


1 


2 X 10" 


5 


0.30 


11.62 


12.42 


13.20 


He32 


1 


2 X 10" 


5 


0.32 


11.32 


12.05 


12.78 


Dl 


1 


1 X 10" 


5 


0.25 


12.71 


13.74 


14.55 


D2 (He25) 


1 


2 X 10" 


5 


0.25 


12.70 


13.55 


14.35 


D4 


1 


4 X 10" 


5 


0.25 


12.45 


13.09 


13.93 


ZO 





2 X 10" 


5 


0.22 


14.37 


15.40 


16.30 


ZO.Ol 


0.01 


2 X 10" 


5 


0.22 


13.59 


14.56 


15.39 


ZO.l 


0.1 


2 X 10" 


5 


0.22 


12.92 


13.82 


14.65 


Z0.3 


0.3 


2 X 10" 


5 


0.22 


12.48 


13.33 


14.14 


Zl.O'^ 


1 


2 X 10" 


5 


0.22 


12.00 


12.79 


13.56 


T0.3 


10-1/2 


1 X 10" 


5 


0.25 


12.92 


13.77 


14.67 


Tl 


1 


2 X 10" 


5 


0.25 


12.54 


13.48 


14.33 


T3 


101/2 


4 X 10" 


5 


0.25 


12.20 


13.13 


13.86 



Models He22-He32 and D1-D4 are calculated with a IBurrows et al.l l| 19971 ) radiative boundary condition Models 
ZO.Ol-Zl.O are calculated with opacities from lAllard fc Hauschildtl (| 19951 '). incorporated in atmospheric boundary 
conditions as described in IBurrows et al] (|2001l ) (D. Saumon, private communication). Models T0.3-T3 have radiative 
boundary conditions calculated with COOLTLUSTY. Model ZO's radiative boundary condition comes from I Saumon et al.l 
(|T99l ). 

^Our adopted value of Zq — 0.0189 is taken from lAnders fc Grevessd (|1989l ). This is in contrast to the 0.0142 value 
from the more recent work of , Asplund et all (j2009l ). In absolute units, the quantities in the table correspond to the 
following values: x = 0; 0.01 x .Zq = 1.89 x 10"*; 0.3 x = 0.00598. 

''This model shows up in two rows, named both He25 and D2, and is one of two "fiducial" models. 
'^This is the other "fiducial" model. 
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TABLE 2 

Approximate Edge-Mass Derivatives 



Model Sequence 


Derivative 


Derivative Value 


He22-He32 


(dM/dY) /lOO 


-0.2Mj 


D1-D4 


(dM/dyD,i)/10^ 


-0.16Mj 


D1-D4 


dM/dlogio[yD,il 


-0.8Mj 


ZO.Ol-Zl.O 


dM/dlog,o[Z] 


~0.9Mj 


T0.3-T3 


dM/dlog,o[Z] 


~0.7Mj 



TABLE 3 

Edge-Masses for "Realistic" Treatment of Z 



Metallicity 


M (10%) 
(Mj) 


M (50%) 
(Mj) 


M (90%) 
(Mj) 


0.5 X solar'' 


12.26 


13.10 


13.93 


1 X solar 


11.89 


12.73 


13.57 


2 X solar 


11.39 


12.23 


13.06 


3 X solar 


10.99 


11.83 


12.67 



We start with Y — 0.2703 and add 0.0142 in Y for every unit in Z. We perform a cubic interpolation in Y and use 
the derivatives recorded in the T0.3-T3 row of Table [2] to approximate the combined affects of metallicity on pressure 
support and on atmospheric opacity. 



^Our adopted value of Zq = 0.0189 is taken from lAnders fc Grevessd (|1989f ). In absolute units, the quantities 
in the table correspond to the following values: 0.5 x Zq ^ 0.00945; 2 x Zq = 0.0378; 3 x Zq = 0.0567. 
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Fig. 1. — Evolutionary trajectories as functions of mass for D2 models. Top left: Radius vs. time. Higher mass objects cool 
faster and have smaller radii by 10 Gyr, even though they begin with larger radii. The 15-A/,j and 20-Mj models cool and shrink 
particularly rapidly after the conclusion of deuterium burning. Top Right: Effective temperature vs. time. Bottom left: Ratio of 
nuclear power to total luminosity vs. time. Bottom right: Fraction of initial deuterium burned vs. time. 
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Fig. 2. — Ratio of nuclear power to total luminosity vs. time (left) and fraction of initial deuterium that burns vs. time {right), 
for 13-Mj models. Note that very little deuterium burning occurs after ~1 Gyr. Models are grouped by varying the helium 
fraction (top), the initial deuterium fraction [middle), and the metallicity {bottom). 



13 




Fig. 3. — Fraction of initial deuterium that burns within 10 Gyr vs. Mass (left); and dependence of deuterium-burning edge on 
Y, yu,i, and Z, for different "edge-criteria" {right). 

Left: For a variety of models, the total fraction of the initial deuterium abundance that combusts through nuclear fusion within 
10 Gyr is shown as a function of the object's mass. Note, per Fig. El that if any appreciable fraction of the initial deuterium 
ends up burning, this happens within the first 1 Gyr. Models He25/D2 and Zl.O are shown in all three panels. In each panel, 

horizontal dashed lines are plotted at 10%, 50%, and 90%. 

Top: For solar metallicity and initial deuterium number fraction of 2 X 10~^, models from [Burrows et al.l 1)19971 ) are shown for 
different helium mass-fractions. Greater helium fraction leads to deuterium burnin g at a lower ma. ss. 

Middle: For solar metallicity and helium abundance (by mass) of 0.25, models froni I Burrows et al.l 119971) are shown for different 
initial deuterium abundances. Greater initial deuterium abundance leads to deuterium burning at lower mass. 
Bottom: For inital deuterium abundance of 2 X 10~^ and helium abundance of 0.25, models using an AUard & Hauschildt { 1993) 
boundary condition (ZO.Ol, ZO.l, Z0.3, Zl.O) are shown. A model with zero metallicity (ZO, Saumon et al. 1994) is also shown. 
Greater metallicity leads to deuterium burning at lower mass. 

Right: In each panel, the red, green, and blue curves correspond to edge-criteria of 10%, 50%, and 90% of the initial deuterium 
burning within 10 Gyr. These correspond to the masses at which the fraction-burned curves in the left panel of this figure cross 
the three horizontal dashed lines. Individual models are coded by the same colors as in the left panel. The mass of the deuterium 
burning edge is shown as a function of helium fraction (top), initial deuterium fraction (middle), and metallicity (bottom). The 
"fiducial model" He25/D2 is represented (in the top two panels) with a large black square, filled with yellow. 



